Microscopic study of ^^°Pu, mean-field and beyond 
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The influence of exact angular-momentum projection and configuration mixing on properties of a 
heavy, well-deformed nucleus is discussed for the example of ^'"'Pu. Starting from a self-consistent 
model using Skyrme interactions, we analyze the resulting modifications of the deformation energy, 
the fission barrier height, the excitation energy of the superdeformed minimum associated with the 
fission isomer, the structure of the lowest rotational bands with normal deformation and superde- 
formation, and the corresponding quadrupole moments and transition moments. We present results 
obtained with the Skyrme interactions SLy4 and SLy6, which have slightly different surface tensions. 

PACS numbers: 21.60.Jz, 21.30.Fe, 21.10.-k, 27.90. -t-b 



I. INTRODUCTION 

Microscopic mean-field methods 0| are particularly 
well suited to describe nuclei with a well defined shape. 
When the energy of a nucleus depends softly on a shape 
degree of freedom or presents several minima as a func- 
tion of this shape, correlations beyond mean field can af- 
fect the properties of the ground state strongly. In such 
cases, the two most relevant types of correlations are as- 
sociated with the rotation of the nucleus and with its 
vibrations with respect to deformation. The inclusion of 
rotational correlations can be performed by a symmetry 
restoration and that of vibrations by a mixing of mean- 
field states corresponding to different shapes. In both 
cases, this requires to go beyond mean-field models. 

We have recently developed a method which achieves 
these goals 0. Applications have been carried out for 
neutron-deficient Pb isotopes UQ. The low-energy spec- 
trum of these nuclei varies rapidly with neutron num- 
ber with states exhibiting strong mixing between oblate, 
spherical and prolate configurations. Qualitative prop- 
erties of their spectra, including transition properties, 
were nicely explained. However, since the results depend 
strongly on the amount of mixing between several con- 
figurations, a detailed agreement with the data has not 
been achieved. 

A very different situation occurs when the mean-field 
approximation is better justified as a first approximation, 
such as when coexisting states lie in well separated en- 
ergy minima. This is the case at low excitation energies 
for superdeformed bands in nuclei around Hg and for fis- 
sion isomers 0. Nevertheless, it is only using beyond 
mean-field models that one can calculate spectra with 
well defined spin assignments as well as the correspond- 
ing transition probabilities. 

The nucleus ^"'"Pu has often been used as a benchmark 
to study mean-field theories and effective interactions. 
We present here an application of our method to this 
nucleus. It will allow us to address the following issues: 

(i) are quadrupole correlations infiuencing a well- 



deformed nucleus a priori well-described by mean-field 
calculations? 

(ii) how does the exact angular momentum projection 
modify the fission barrier and the excitation energy of 
fission isomers? 

(iii) how much do the predictions depend upon differ- 
ent parameterizations of the effective interaction? 

In what follows, we shortly recall the basic ingredients 
of the theory, then we present our results for the spectrum 
of ^^"Pu in the ground state and the superdeformed well. 
The fission barrier obtained with two different effective 
interactions is discussed and compared with earlier, more 
phenomenological, approaches. 



II. THE MODEL 

The starting point of our method is a set of HF-I-BCS 
wave functions generated by self-consistent mean-field 
calculations with a constraint on a collective coordi- 
nate, the axial quadrupole moment q — {Q20) in the 
present study. In the language of the spherical nuclear 
shell model, such mean-field states incorporate particle- 
particle (pairing) correlations as well as many-particle 
many-hole correlations by allowing deformations of the 
nucleus in its intrinsic frame. As a consequence, the 
mean-field states break several symmetries of the exact 
many-body states. This symmetry violation makes it dif- 
ficult to relate mean-field results to spectroscopic data 
which are obtained in the laboratory frame of reference. 
The second step of our method is a restoration of the sym- 
metries associated with particle numbers and rotation. 
Another ambiguity in the interpretation of mean-field 
results arises from the non-orthogonality of mean-field 
states corresponding to different quadrupole moments, 
so that different minima in a potential landscape cannot 
always be safely associated with different physical states. 
This difficulty is resolved in the third step of our method 
by variational mixing of symmetry-restored mean-field 
states corresponding to different quadrupole moments. 
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The method that we use is a discretized version of the 
generator coordinate method. It removes the contribu- 
tion of vibrational excitations to the ground state and, 
at the same time, permits to construct a spectrum of 
excited states. 

In our method, the same effective interaction is used to 
generate the mean-field states and to perform the config- 
uration mixing. We present below results obtained with 
two different Skyrme interactions, SLy4 and SLy6 In 
both cases a density-dependent zero-range interaction is 
used in the pairing channel. We use the same strength 
as in previous studies, —1250 MeV fm'^ and two cutoffs, 
above and below the Fermi energy, as defined in Ref. . 
The two Skyrme parameterizations differ mostly by their 
surface tension: the surface energy coefficient obtained 
from Hartree-Fock calculations of semi-infinite nuclear 
matter is lower for SLy6, 17.74 MeV, than for SLy4, 
18.37 MeV j^j- Such a difference is expected to affect 
significantly the deformation energy at large quadrupole 
moments. Both parameterizations have been fitted in an 
identical way. Their differences have their origin in a 
different choice for the treatment of the spurious center- 
of-mass motion (cm.): a fully variational cm. one for 
SLy6 and a simpler one-body approximation for SLy4. 
The energy differences due to these two schemes induce 
slight differences in the properties of the interactions, see 
Ref. [9| for a detailed discussion. 

As our main goal is an investigation of the overall ef- 
fect of symmetry restoration and configuration mixing in 
a heavy, well-deformed nucleus, we restrict ourselves to 
axial and reflection-symmetric shape degrees of freedom. 
It has been shown that the fission barrier height obtained 
with SLy6 is in agreement with experiment within 1 MeV 
when octupole deformation is taken into account fll|. 

Our method has many interesting properties. Its 
sole phenomenological ingredient is the effective nucleon- 
nucleon interaction, which has been adjusted once and 
for all on generic nuclear properties. From a numerical 
point of view, it is simple enough to be applied through- 
out the mass table up to superheavy nuclei, utilizing the 
full model space of single-particle states with the proper 
coupling to the continuum. Another attractive aspect of 
the method is that it allows to determine electric tran- 
sition probabilities directly in the laboratory frame be- 
tween any pair of states. Finally, the method has the 
advantage that its results can be interpreted within the 
intuitive picture of intrinsic shapes and shells of single- 
particle states that is offered by the framework of mean- 
field models. More details on the method can be found 
in Refs. 

Spectroscopic quadrupole moments and B{E\) values 
are determined directly in the laboratory frame of refer- 
ence 0. To connect our results with other approaches, 
it is interesting to derive quantities analogous to intrin- 
sic frame parameters from spectroscopic or transition mo- 
ments. An intrinsic charge quadrupole moment (5^2 {J^ ^) 
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FIG. 1: Deformation energy curve of ^''"Pu obtained with 
SLy4 projected on A'^ and Z (dashed line) and projected on 
A^, Z and J = (solid line). All energies are normalized to 
the deformed ground-state value of each curve. The available 
experimental data for the excitation energy of the superde- 
formed band head are shown at arbitrary deformation (see 
text). Shapes along the path are indicated by the density 
contours at p = 0.07 fm~^. 



can be determined from B{E2) values: 



'I67r B{E2,J ^ J -2) 
~(J020|J-20)2e2 



(1) 



or can be related to the spectroscopic quadrupole mo- 
ment Qc{J, k) via the relation 



2J + 3 



,{J,k). 



(2) 



We also adopt the sharp edge liquid drop relation 
to relate the (32 deformation parameter and the axial 
quadrupole moment Q2 



P2 



5 AttQ2 
167r 3i?2A ' 



(3) 



where the nuclear radius R in fm at zero deformation is 
related to the mass A according to the standard formula 
R = 1.2^1/3. 



III. RESULTS 

A. Deformation energy 

There is a large set of data on fission barriers of ac- 
tinide nuclei ,11. 12J . Among them, the double- humped 
fission barrier of ^'"'Pu has been used as a benchmark 
for mean-field models and effective interactions. First 
calculations were performed with Skyrme forces and the 
Hartree-Fock-|-BCS method with the Gogny force 
and the HFB method 0| or with relativistic Lagrangians 
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FIG. 2: Same as Figure for the interaction SLy6. 

and the relativistic mean- field method (RMF) [H. Semi- 
classical approximations of the mean-field method were 
performed with different Skyrme forces in ref Il7j |. 
Several Skyrme interactions and RMF Lagrangians were 
compared to the data in [Ts^ . Axial and triaxial barriers 
obtained with Skyrme, Gogny and RMF forces are com- 
pared in Ref. |l| . Finally, the excitation energy of fission 
isomers has been studied with a variety of Skyrme forces 
in Ref. 

The deformation energy curves obtained after particle- 
number projection and particle-number -I- angular- 
momentum projection on J = are presented in Figures 
Hand El for SLy4 and SLy6 interactions. For all curves, 
the energy of the ground state is taken as zero. The 
ground state and the fission isomer after projection are 
obtained from the mean-field minima. With the normal- 
ization that we have chosen to plot the results, the gain 
of energy obtained for the ground state by angular mo- 
mentum projection is given by the difference between the 
curves at spherical shape. This gain is around 3.0 MeV 
for both interactions, bringing the calculated total energy 
closer to the experimental one. 

The fission isomer is obtained at a (32 value around 0.9; 
angular-momentum projection lowers its excitation en- 
ergy by about 1 MeV for both forces, from 4.3 to 3.3 MeV 
for SLy4, and from 2.6 to 1.6 MeV for SLy6. 



B. Rotational energy 

Angular momentum projection provides the exact cor- 
rection for the spurious rotational energy of the mean- 
field states. It is given by the difference between the 
binding energies before and after projection on angular 
momentum J = 



(4) 



This difference does not depend much on the Skyrme pa- 
rameterization and is plotted in the lower panel of Fig. El 
for SLy4. It is zero at spherical shape and increases first 
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FIG. 3: Top to bottom: the average angular momentum I of 
the mean-field states obtained from ( J'^) = fi? the mo- 

ment of inertia calculated from the difference of the projected 
J — 2 and J — energy curves (solid line) and the Belyaev 
moment of inertia (dotted line); and the rotational energy 
obtained from the energy difference between the mean-field 
and the J = energy curves (solid line) and the rotational 
correction, Eq. @. 



rapidly to values around 3 MeV for deformations smaller 
than 1/32 1 < 0.1, and then moderately for larger deforma- 
tions. A similar behavior has been obtained in most of 
our previous calculations II, [13. The topology of the fis- 
sion barrier is not much affected by angular momentum 
projection. The height of the second barrier is decreased 
by about 800 keV with respect to the fission isomer and 
by 1.5 MeV with respect to the ground state. 

Rotational corrections are sometimes incorporated 
phenomenologically as a perturbation to mean-field cal- 
culations. In particular, it is an ingredient in the mass 
formulae based on Skyrme forces 20]. In the same way, 
the fission barrier of ^''"Pu including a rotational cor- 
rection was used as a constraint in the fit of the Gogny 
interaction jljj. In both cases, the rotational energy has 
the form 



2e(/32)' 



(5) 
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FIG. 4: Collective wave functions (upper panel) of the lowest 0"*", 2^, and 4"'" states. The lower panels display the corresponding 
excitation energies at the average deformation of the mean-field states from which they are built, together with the projected 
energy curve. 



where {J"^) is the mean value of the square of the angular 
momentum for the mean-field state, and the moment of 
inertia O is determined from an approximate cranking 
formula 

6Bclyacv = 2 ^ ~ "^'"j)- 

The sum in Eq. ^ runs over the single-particle states 
\i) and \j) in a given deformed mean-field state, and Ei 
and Ej are the corresponding quasi-particle energies. For 
the large-scale mass fits of Ref. (201 , the actual moment 
of inertia taken is a mixture of Eq. © and of the rigid- 
body moment of inertia. More involved approximations 
for the moment of inertia have been developed, see e.g. 
[21L and references therein, but are rarely used. We 
can extract a moment of inertia from our projected mean- 
field calculations using the energy difference between the 
energy curves for J = and J = 2 as 

In Fig. |3| we compare the exact and approximate rota- 
tional energies as a function of deformation. We also 
show the angular momentum dispersion of the mean- 
field wave functions and the moment of inertia given by 
Eqn. © and obtained from the 2+ excitation energy. 
Both of these moments of inertia are rather close for de- 
formations around that of the fission isomer. For lower 
deformations, the Belyaev moment of inertia is signifi- 
cantly larger than the "exact" one. For larger deforma- 
tions, the Belyaev moment does not increase as rapidly 
as the one extracted from our calculation. The disper- 
sion of the angular momentum of the mean-field wave 
function increases also with deformation in a way which 



is partly compensated by the increase of the Belyaev 
moment of inertia. However, this compensation is not 
strong enough and the Belyaev rotational energy cor- 
rection overestimates at large deformations the energy 
correction obtained by exact projection. Note that this 
correction varies only by 1.5 MeV from the deformation 
corresponding to the ground state to the external fission 
barrier. 



C. Configuration mixing 

The properties of the four lowest states obtained from 
the configuration mixing calculation are given in Tabled 
for SLy4 and Table ITU for SLy6. The corresponding col- 
lective wave functions are shown in Figure^lfor the SLy4 
interaction. The states separate nicely into four rota- 
tional bands, two located in the prolate normal-deformed 
(ND) minimum, and two in the superdeformed (SD) one. 
The wave functions are all confined within either the ND 
or the SD wells. States obtained with SLy6 are similar, 
except for an overall shift in the excitation energies of 
the states located in the superdeformed well; the SD 0+ 
band-head has an excitation energy of 2.99 MeV with 
SLy4 and 1.25 MeV with SLy6. The actual experimen- 
tal value is in between the two, although there are some 
conflicting values given in the literature 2.4 ± 0.3 MeV 
[HI, « 2.8 MeV [3 and 2.25 ± 0.20 MeV 

The band head of the second ND band is obtained 
at 3.79 MeV for both forces. Similarly, the excitation 
energy of the second SD band with respect to the first 
SD band head is 1.39 MeV for SLy4 and 1.27 MeV for 
SLy6. This suggests that the excitation energies within a 
well are fairly independent from the surface tension of the 
Skyrme force. Both interactions give also similar values 
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state 


E Qs Qo"^ 
(MeV) (eb) (eb) 




B{E2)l 


(eb) 


4** 


0^ 
< 


0.000 
0.083 
0.277 


-3.4 
-4.3 


11.9 
11.9 


0.300 
0.300 


2.80 
4.00 


11.9 
11.9 


0.300 
0.300 


O3+ 
23+ 
43+ 


3.793 
3.880 
4.088 


-3.4 
-4.3 


11.9 
12.0 


0.301 
0.303 


2.82 
4.07 


11.9 
12.0 


0.301 
0.302 




2.953 
2.978 
3.045 


-10.3 
-13.1 


36.0 
36.0 


0.911 
0.911 


25.8 
36.8 


36.0 
36.0 


0.911 
0.911 


24+ 
44+ 


4.338 
4.364 
4.429 


-10.4 
-13.3 


36.5 
36.5 


0.922 
0.922 


26.4 
37.8 


36.5 
36.5 


0.922 
0.922 



TABLE I: Properties of the rotational bands of ^"'"Pu 
obtained with SLy4: excitation energy E, spectroscopic 
quadrupole moment Qs, corresponding quadrupole moment 
Qq'^ and dimensionless deformation /Sj"' in the intrinsic frame, 
reduced E2 transition probability B{E2)l, and corresponding 

quadrupole moment Qq*' and dimensionless deformation /Jj*' 
in the intrinsic frame. 



state 


E Qs Qo^' 
(MeV) (e b) (e b) 




B{E2)l 
(e^b^) (eb) 




ot 

2^ 
4t 


0.000 
0.083 
0.273 


-3.4 
-4.3 


11.9 
11.9 


0.300 
0.301 


2.81 
4.02 


11.9 
11.9 


0.300 
0.301 


Ot 
2t 
44+ 


3.794 
3.879 
4.082 


-3.4 
-4.4 


12.0 
12.1 


0.304 
0.306 


2.88 
4.15 


12.0 
12.1 


0.304 
0.305 


Ot 
2J 


1.251 
1.277 
1.338 


-10.4 
-13.2 


36.3 
36.3 


0.916 
0.917 


26.1 
37.4 


36.2 
36.3 


0.915 
0.916 


ot 

2t 

4 


2.519 
2.550 
2.611 


-10.5 
-13.4 


36.7 
36.7 


0.928 
0.928 


27.0 
38.3 


36.8 
36.7 


0.931 
0.928 



TABLE IL The same as Table ID but for SLy6. 



for the excitation energies within the bands. 

The gro und state band is known up to very high spin 
|2^|2a,|23. It has been suggested that static octupole de- 
formation plays a role to explain the behaviour at large 
angular momentum Ref . ,2^ . States below the 6+ 
decay mainly by internal electron conversion, so the cor- 
responding transitions have not been detected in 7-ray 
spectroscopy. 

The lowest levels in the ground-state band are reported 
in the NUDAT data base ^ at 42.824 (2+), 141.690 
(4+), and 294.319 keV (6+). Our calculation overesti- 
mates these energies by almost a factor two. The exper- 
imental energies in the SD well are 20.1 keV for the 2+ 



and 66.8 keV for the 4+ for the SDl band, and 769.9 
(0+), 785.1 (2+), and 825.0 keV (4+) for the SD2 band. 
They are also overestimated by our model. On the con- 
trary, the calculated B{E2; 2+ 0+) value is in excel- 
lent agreement with the experimental one of 26660 ± 360 

fm"' obtained from Coulomb excitation 30j. 

A similar overestimation of excitation energies has 
been found for other nuclei with our model ^] and 
in a similar framework using a Gogny interaction |3l| . 
A hint on a possible origin of this discrepancy can be 
found from cranked mean- field calculations. In this case, 
neither configuration mixing nor restoration of symme- 
tries are performed. However, time- reversal invariance 
is broken and the mean-field potential is optimised for 
each J value and not only for J = 0. An unprojected 
cranked HFB calculation i^S] using the same effective in- 
teraction as here gives excitation energies for the ground- 
state band of 0.030 (2+), 0.121 (4+) and 0.271 MeV (6+), 
respectively, in much better agreement with the data. 
Looking at Figure^ one can see that, within a band, the 
amplitudes of the collective wave functions for different 
J values differ by less than 1% for each deformation. The 
wave functions for different J values are therefore prob- 
ably too close to each other in our model. The slight 
breaking of time-reversal symmetry of mean-field states 
subject to a cranking constraint might be sufficient to 
provide a better starting point for exact projection and 
configuration mixing for the J different from states. 

The deformation of the ground state stays remarkably 
constant at all levels of approximations: from (32 — 0.29 
for the minima of the mean-field and projected energy 
curves to /3^'-'(2+) = 0.30 as deduced from Eqn. lEJl for 
both the spectroscopic and the transition quadrupole mo- 
ments. All these deformations agree with the one de- 
duced from the experimental B{E2) value, f32 = 0.29. 

Since we obtain nearly equal /^j"'' and ^ values all 
along the bands, the use of the rotor model is well justi- 
fied to describe the four bands. 

We obtain also very large EO transition matrix ele- 
ments between states in the same well. With SLy4, 
the M{EO;Ot Ot) in the first well has a value of 
-29 e fm^ corresponding to p{EO) = —0.52, while the 
the M{EO;Ot Ot) in the second weh is 87 e fm^, 
or p{EO) = 1.6. As our model predicts the collective 
wave functions of all members of a rotational band to be 
very similar, the transition moments are calculated to be 
very similar as well for all EO transitions within a band. 
The values obtained with the SLy6 interation differ only 
marginally from those of SLy4. The -EO matrix elements 
that we obtain for transitions between states in the SD 
and the normal-deformed well are very small. 

IV. DISCUSSION AND SUMMARY 

The inclusion of correlations beyond mean-field con- 
firms that many properties of the ^"^^Pu nucleus are al- 
ready well described at the mean-field level of approxima- 
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tion. The overall structure of the potential energy curve 
is not altered by angular momentum projection, with a 
well defined prolate ground state and a fission isomer in 
narrow potential wells. Thanks to that, the properties of 
the lowest state in each well after configuration mixing 
are very close to the properties of the mean-field min- 
ima. The ground-state and fission isomer wave functions 
have a Gaussian shape and do not spread much around 
their respective minima. The superdeformed minimum 
in the potential energy surface is, however, too wide to 
confine completely the wave functions of the lowest SD 
band which are more spread, but still of Gaussian shape. 

As expected, the total binding energy is slightly in- 
creased by angular momentum projection, and thereby 
comes closer to the experimental value. Compared to 
the ground state, angular-momentum projection lowers 
the (axial) inner barrier by about 0.6 MeV, the fission 
isomer by about 1 MeV, and the (reflection-symmetric) 
outer barrier by about 2 MeV. These changes of the 
potential landscape are going in the same direction as 
a decrease of the surface tension of the effective mean- 
field interaction and should therefore be considered when 
predicting superdeformed band heads and fission barrier 
heights. The schematic rotational correction used in the 
literature gives a too large reduction of the fission barrier 
heights by at least 2 MeV. 

The energy of the superdeformed fission isomer turns 
out to be too low with SLy6, while for SLy4 it is slightly 
too high. We obtained similar results for SD band heads 
of Pb isotopes in the A w 190 region 0. This suggests, 
that the surface tension might be too low for SLy6, while 



it is slightly too high for SLy4. This is apparently in con- 
tradiction with Ref . i^J , where it was argued on the basis 
of pure mean-field calculations for ^^°Pu, that the surface 
tension of SLy6 is more realistic than that of SLy4. How- 
ever, surface tension is not the only ingredient responsi- 
ble for the energy of superdeformed states, as their exis- 
tence is usually caused by a shell effect. Unfortunately, it 
is hard to disentangle the contribution from the macro- 
scopic properties of the forces from the shell structure 
which is also not identical for SLy4 and SLy6 at large 
deformations. 

The description of excited states in both wells is only 
partly satisfactory. While we obtain quite good results 
for E2 transition moments and deformations, neither the 
excitation energies of the excited 0+ band heads nor the 
excitation energies within the bands are well reproduced. 
Unprojected cranked HFB calculations without any addi- 
tional correlations perform much better for in-band tran- 
sitions. This strongly suggests to extend our method to 
the use of such states as a starting set of wave func- 
tions. Such a development is currently underway. Using 
cranked mean-field states, however, will not change the 
too large excitation energies of 0+ band heads. 
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